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Abstract 

In return for the long-standing contributions of Physics to Biology, now the inverse way is 
frequently traveled through in order to think about many physics phenomena. In this vein, evo- 
lutionary algorithms, particularly genetic algorithms, are being more and more used as a tool to 
deal with several Physics problems. Here, we show how to apply a genetic algorithm to describe 
the homogeneous electron gas. 
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Several methods, techniques and ideas taken from Physics have been of primary impor- 
tance for the notable development of Biology in the last sixty years. To agree with this 
statement is it enough to think about the application of X Ray crystallography to the reso- 
lution of protein tertiary and quaternary structures, or the fundamental role that Quantum 



lave played in the establishment of 
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However, since some years ago, 



Mechanics in general and Molecular Physics in particular 
modern (Molecular) Biology, as predicted by Schrodinger 
the inverse way is also frequently traveled by physicists. Biological methods and concepts, 
as well as its terminology itself, are more and more used in Physics. 

In this context, in the last decades, there has been an increasing interest in algo- 
rithms based on the Darwinian evolution principle!^. Genetic algorithms evolution- 
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ary programming p|,p|, game-playing strategies |2I| and genetic programming have found 
a wide field of applications, not just in Physics but also in other areas where optimization 
plays an important role, such as financial markets, artificial intelligence, etc. In particular, 
genetic algorithms tackle even complex problems with surprising efficiency and robustness. 
In Physics they have been used in calculations thet involve from single Schrodinger particles 
in diverse potentials to astrophysical systems, running through lattice systems, spin glass 
models, molecules and clusters. The differences among the various evolutionary algorithms 
can be found not in the basic principles but in the details of the selection, reproduction and 
mutation procedures. 

In general, an evolutionary algorithm is based on three main statements: 

a) It is a process that works at the chromosomic level. Each individual is codified as a 
set of chromosomes. 

b) The process follows the Darwinian theory of evolution, say, the survival and reproduc- 
tion of the fittest in a changing environment. 

c) The evolutionary process takes place at the reproduction stage. It is in this stage when 
mutation and crossover occurs. As a result, the progeny chromosomes can differ from their 
parents ones. 

Starting from a guess initial population, an evolutionary algorithm basically generates 
consecutive generations (offprints). These are formed by a set of chromosomes, or character 
(genes) chains, which represent possible solutions to the problem under consideration. At 
each algorithm step, a fitness function is applied to the whole set of chromosomes of the 
corresponding generation in order to check the goodness of the codified solution. Then, 
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according to their fitting capacity, couples of chromosomes, to which the crossover operator 
will be applied, are chosen. Also, at each step, a mutation operator is applied to a number 
of randomly chosen chromosomes. 

The two most commonly used methods to randomly select the chromosomes are: 

i) The roulette wheel algorithm. It consists in building a roulette, so that to each chro- 
mosome corresponds a circular sector proportional to its fitness. 

ii) The tournament method. After shuffling the population, their chromosomes are made 
to compete among them in groups of a given size (generally in pairs). The winners will be 
those chromosomes with highest fitness. If we consider a binary tournament, say the com- 
petition is between pairs, the population must be shuffled twice. This technique guarantees 
copies of the best individual among the parents of the next generation. 

After this selection, we proceed with the sexual reproduction or crossing of the chosen 
individuals. In this stage, the survivors exchange chromosomic material and the result- 
ing chromosomes will codify the individuals of the next generation. The forms of sexual 
reproduction most commonly used are: 

i) With one crossing point. This point is randomly chosen on the chain length, and all 
the chain portion between the crossing point and the chain end is exchanged. 

ii) With two crossing points. The portion to be exchanged is in between two randomly 
chosen points. 

For the algorithm implementation, the crossover normally has an assigned percentage that 
determines the frequency of its occurrence. This means that not all of the chromosomes will 
exchange material but some of them will pass intact to the next generation. As a matter 
of fact, there is a technique, named elitism, in which the fittest individual along several 
generations does not cross with any of the other ones and keeps intact until an individual 
fitter than itself appears. 

Besides the selection and crossover, there is another operation, mutation, that produces a 
change in one of the characters or genes of a randomly chosen chromosome. This operation 
allows to introduce new chromosomic material into the population. As for the crossover, 
the mutation is handled as a percentage that determines its occurrence frequency. This 
percentage is, generally, not greater than 5%, quite below the crossover percentage. 

Once the selected chromosomes have been crossed and muted, we need some substitution 
method. Namely, we must choose, among those individuals, which ones will be substituted 
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for the new progeny. Two main substitution ways are usually considered. In one of them, all 
modified parents are substituted for the generated new individuals. In this way an individual 
does never coexist with its parents. In the other one, only the worse fitted individuals of the 
whole population are substituted, thus allowing the coexistence among parents and progeny. 

Since the answer to our problem is almost always unknown, we must establish some 
criterion to stop the algorithm. We can mention two such criteria: i) the algorithm is run 
along a maximum number of generations; ii) the algorithm is ended when the population 
stabilization has been reached, i.e. when all, or most of, the individuals have the same 
fitness. 

In this letter we apply some of the previous ideas in order to calculate the pair distribution 
function (PDF) of an homogeneous electron gas. For clarity we consider the one-dimensional 
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version of the electron gas, but the method is directly extended to higher dimensions 

n 

Let us consider the Sommerfeld-Pauli model of metallic solids[li| in ID, namely N elec- 
trons of mass m moving along a segment of the axis x of length L . The position of the ith 
electron is denoted with Xi {i = 1,2, ...,N). Electrons i and j interact one with the other 
through a pair potential v (xij) pair potential is often used to model 

a quantum wire of width 6. The system Hamiltonian then reads: 
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For the A^-body wave function we use the trial form 



(1) 



^ {x) = exp 

i<j i<j 

where ~x = {xi,X2, ■■■,X]\f) denotes the system configuration and 
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271 p J So{k) 
Here p = N/ L is the system number density and go {x) is the ideal {v (x) = 0) electron gas 
pair distribution function with its Fourier transform 5*0 (k), the structure factor. 

To apply the algorithm, we take M distances yi = 0, 1/2 = d,...,yM = {M — l)d so 
that, for any the function u{xij) is represented by the random string u{xij =yi), 
u {xij = y2),..-,u {xij = yu), where each u {xij = yu) is a random real number 7^^- G [0, 1] 



(rounded to an established number n of decimals): u{xij = yk) = 'Jij. Our initial popula- 
tion is then formed by Np random replicas of each of the (A^ — 1) /2 strings. In this work 
we have taken = 100, M = 99, = 3 and Np = 100. 

Given a string 'jljy'Jij, ■■■,'Jij , the encoding consists in replacing the sequence of real 
numbers G [0, 1] by a single natural number obtained by putting their decimal parts 
one next to the other. For example, taking n = 3: the string 0.137, 0.935, 0.466 gives 
137935. ..466. In genetic terms, the encoding produces the chromosomic structures of the 
parents (strings). The inverse process is called decoding. 

For the parents population, the energy of the ath replica (a = 1, 2, ...Np) is given by 

where for the mean kinetic energy we use the Jackson- Feenberg formulae |l2j|. 

In order to calculate the fitness of the replica a, we previously need to estimate the energy 
error: ta = \{H — E^) '^a\ and define the fitness function as fa = e""". A solution is reached 
when /a ^ 1. 

The calculation proceeds by dividing the population of Np replicas into A'p/2 couples. 
The couples are randomly chosen by using the roulette wheel algorithm [? ]. This is done 
by defining the sums F = Yl!l=ifa and = J2a=i W = l^^, ...,Np). Then, a random 
number r G [0, F] is generated and the unique index 6 such that Ss^i < r < Ss is picked up. 

Once the first generation of replicas (parents) has been generated and divided into cou- 
ples, the second generation (offspring) can be generated by applying the crossover operator 
between the members of each one. Sometimes, some of the members of the new replicas 
generation can be changed by applying the mutation operator. 

Given a couple of replicas, the crossover operator is defined by generating a new random 
number c G [0, 1] which is compared with a pre-established crossover probability p G [0, 1]. 
If c < p, the crossover operator acts by interchanging all the digits from the sth position to 
the end of the replica between the members of the couple. Here s is a random integer such 
that 1 < s < nM. for example if the couple is 

153280472. ..337 

768325399. ..069 

and s = 4, the new offspring couple will be 
153225399. ..069 
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768380472... 337. 

To apply the mutation operator we first randomly select those offsprings that will mutate. 
Then, for each of these offsprings, a gene (a digit) randomly chosen is changed by a random 
integer number i G [0, 9] . 

The algorithm was stopped after 2 x 10^ generations. Then, the PDF at the distances 
yi,y2, ■■■,yM is calculated 



, , N{N — 1) J dxsdxi. . .dxN\'^{xi,X2, . . .xn]\ 

g yxj = 2- (, / 

p J dxidx2 . ■ . dxN \'^{xi,X2, ■ ■ ■ Xn)\ 

The so-calculated PDF's were compared with those obtained from Fermi hypernetted 
chain (FHNC) and Monte Carlo (MC) variational approaches with good agreement for dif- 
ferent values of the relevant parameters (Figs. 1-3). FHNC and MC probably are the two 
most reliable methods available to calculate the electron gas correlations. 

Taking into account the good agreement among the curves, and the above mentioned 
efficiency and robustness of the genetic algorithm, we conclude that the application of this 
kind of algoHth,„s to many-body problems dese.yes some attention from physicists m order 
to achieve a better understanding of the involved physical subtletiesj 
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Figure Captions 



Figure 1. The pair distribution function for the ID homogenous electron gas with 
Wigner-Seitz radius Vs—l and screening parameter 5—1. 
Figure 2. Same as Fig. 1 for rs=l and (5=2. 
Figure 3. Same as Fig. 1 for r5=4.5 and 6=1. 
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Fig.1 Stoico etal. 




Fig. 2 Stoico etal. 




Fig. 3 Stoico etal. 



